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We propose that negative absolute temperatures in ultracold atomic clouds in optical lattices can 
be used to simulate quantum systems in new regions of phase diagrams. First we discuss how the 
attractive SU(3) Hubbard model in three dimensions can be realized using repulsively interacting 
173 Yb atoms, then we consider how an antiferromagnetic S=l spin chain could be simulated using 
spinor 87 Rb or 23 Na atoms. The general idea to achieve negative absolute temperatures is to reverse 
the sign of the external harmonic potential. Energy conservation in a deep optical lattice imposes 
a constraint on the dynamics of the cloud, which can relax toward a T < state. As the process is 
strongly non-adiabatic, we estimate the change of the entropy. 
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The aim of dedicated quantum simulations is to study 
experimentally the properties of a complex quantum sys- 
tem by another quantum system which is less compli- 
cated or, at least, better understood. Among these ap- 
proaches one has to mention graphene, where electrons 
can be described by a 2 + 1-dimensional Dirac equa- 
tion bilayer 3 He, which shows quantum critical be- 
havior Q that resembles heavy fermion materials Q; 
atom corals Q where quantum billiards could be simu- 
lated; or cavity quantum electrodynamics where one can 
study the coupling between light and matter Q . 

Ultracold atomic clouds nevertheless have a distin- 
guished place among quantum simulations. During the 
last two decades, experimental control developed to a 
high level in the tuning of microscopic parameters and in 
detection. Ultracold atomic clouds are also very well iso- 
lated and free of unwanted perturbations, like disorder, 
phonons, etc. Finally, the relevant time and length scales 
are typically larger than in solid-state systems; therefore 
the behavior is more tractable experimentally, even in 
out-of-equilibrium situations. 

A number of experiments with ultracold atomic clouds 
benefited from the advantages discussed above. In ad- 
dition to the realization of the Hubbard model for two 
fermionic species [f| 0] (highly relevant for high transi- 
tion temperature superconductivity Q), ultracold atomic 
clouds have been used to study the Mott transition for 
bosons Q and the Anderson localization by a disordered 
potential [1C[, b ut also more abstract concepts such as 
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gauge fields [llj, [12J and black holes 

In this work, we shall investigate how negative absolute 
temperatures, T < 0, in ultracold atomic experiments 
can grant access to unexplored regions of phase diagrams. 
After a general discussion, we will focus on two model 
Hamiltonians. 

The idea of negative absolute temperatures is neither 
complicated nor new [TJ, [l5| . The reason why T < is 
only seldom discussed is that for most physical systems, 
negative T is not possible in equilibrium. However, as 
proposed in Refs. [H, ultracold atomic clouds in op- 
tical lattices can exhibit T < 0, presenting new possibil- 



ities to study correlated quantum systems. 

To gain some insight about the T < regime, and to 
show connections to the more familiar case of T > 0, 
let us consider the equilibrium partition function of a 
quantum system, 



Z = Tr cxp(-/3H) 



(1) 



Here j3 = 1/T is the inverse temperature (we set the 
Boltzmann constant ks = 1) and H is the Hamiltonian 
with n and E n denoting (many-body) eigenstates and 
the corresponding energies. There is no principle which 
forbids (3 < in general. Nevertheless, most physical 
systems have a spectrum E n that is bounded from below 
but has no upper bound. The simplest example of such a 
spectrum is the kinetic energy of free particles, ~ p 2 , or 
for ultracold atoms the trapping potential m Vqt 2 with 
confinement Vq > 0. It is clear from Eq. ([T]) that without 
an upper bound in energy, no equilibrium is possible at 
/? < as the partition function diverges. Although this 
is the case for the vast majority of physical systems, it 
is not impossible to construct Hamiltonians which have 
only an upper bound in energy, where, in turn, only /3 < 
is possible in equilibrium. We shall below focus on how 
this can be achieved with ultracold atoms. 

From Eq. (JTJ) we can also see that while for /? > the 
lower energy states have higher occupation probabilities 



P(E) 



-I3E 



the opposite, an "inverse population" , is 



realized at /3 < 0. Inverse population of the energy levels 
leads to an unusual condition in matter, which is most 
widely used in lasers [lH ■ 

The most important consequence of Eq. ([T]) is that the 
phase diagram of H with /3 < can be mapped to the 
phase diagram of a system with f3 — —j3 described by 
H = —H. This is the equivalence that could be used 
in quantum simulations to chart new parts of phase di- 
agrams experimentally, as we shall discuss later in the 
context of two models. 

Although establishing the correspondence between the 
Hamiltonians H and — H at ft > and j3 < is com- 
pletely trivial theoretically, it represents a major chal- 



2 



step 1) 



step 2) 



load lattice increase 

lattice depth 




Vg< 0, /3f< 

FIG. 1: (Color online) Reversing the external harmonic po- 
tential, to achieve T < with an ultracold atomic cloud in a 
deep optical lattice. 



lenge in an experiment since ultracold atomic clouds are 
always prepared at T > 0. In order to achieve negative 
T, one has to follow a certain protocol, as discussed in 

Refs. mini, 

which we review here in the case of fcrmions 
in a deep optical lattice. 

From the theoretical aspect, the term which deter- 
mines the sign of (3 is the external potential character- 
ized by Vq, since the kinetic and interaction terms (de- 
scribed by the Hubbard model; see the review by Bloch 
et al. UM) are bounded. Therefore the main step is to 
reverse the external potential Vo — > — Vo and then let the 
cloud relax. The reason why the cloud cannot explode is 
energy conservation: the potential energy of the atoms 
cannot be reduced indefinitely as it can only be converted 
to kinetic and interaction energy of the atoms. Further- 
more, since the external potential is not bounded from 
below, equilibrium at positive temperatures is excluded. 

From the experimental side, the protocol involves the 
following steps (see also Fig. [TJ): (1) Load the cloud to an 
optical lattice in an (insulating) state at Vo > and (3 > 
0; (2) increase the height of the optical lattice suddenly 
to "freeze" the density distributions; (3) reverse the sign 
of the external potential Vo — > — Vq (and possibly adjust 
other control parameters, as well); (4) suddenly reduce 
the height of the optical lattice; and (5) wait for the 
cloud to relax. According to the Boltzmann simulation in 
Ref. [ljj , the cloud should be close to an equilibrium state 
with T < within the typical experimental time scale. 
The protocol works optimally when the initial state is 
an insulator, somewhat restricting the final parameters. 
To access a more extended region of the phase diagram, 
one should change the couplings quasi-adiabatically after 
step 5. 

Although the dynamics restricted by the principle of 
energy conservation can lead to T < in a reversed 
parabolic potential, the final state could be too hot: the 
interesting correlated phases typically found in the low- 
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FIG. 2: (Color online) Schematic phase diagram of the ho- 
mogeneous attractive SU(3) Hubbard model. The quantum 
phase transition between the color superfluid and the trionic 
phase happens when the interaction \U\ is of the order of the 
bandwidth. 



temperature regime might not be accessed experimen- 
tally. The main source of heating is that the protocol 
is inescapably non-adiabatic. Therefore it is important 
to discuss which are the suitable conditions and estimate 
how much entropy is generated. The optimal parameters 
depend on the details of the specific system. Therefore 
we focus on two, rather different but theoretically impor- 
tant models to show how simulations could benefit from 
applying T < 0. 



THE ATTRACTIVE SU(3) HUBBARD 
MODEL 



Fermions hopping between nearest-neighbor lattice 
sites with three possible internal states (colors) and in- 
teracting locally with the same interaction strength are 
described by the SU(3) Hubbard model: 

H su(3) = _ tp [c\ a c ja + H.c] + U J2 »W*tf> (2) 

{ij)a ia<{3 

where tp > is the nearest- neighbor hopping amplitude, 
U is the strength of the local interaction, and c| is a 
fcrmionic creation operator at site i for the fermion color 
a = 1,2,3. 

The model in Eq. with an attractive interaction, 
U < 0, received considerable attention [20t - L23j : This 
Hamiltonian, in addition to being a generalization of the 
two-component case with SU(2) symmetry, also shows 
similarities to quantum chromodynamics, the gauge the- 
ory of quarks describing the strong interaction. In par- 
ticular, both models exhibit a quantum phase transition 
between a superfluid phase with broken color SU(3) sym- 
metry (color superfluid) and a phase where three particles 
form color-singlet bound states (baryons or trions) . If we 
could simulate the attractive SU(3) Hubbard model, we 
could glimpse some puzzles of this quantum field theory 
[24j . A schematic phase diagram is shown in Fig. [5] 

To realize the model in Eq. ^ , the attention was first 
focused on 6 Li atoms in a strong magnetic field, due 
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to the hyperfine structure and the large neg ative back- 
ground scattering length a,b g « — 100nm[25j. Unfortu- 
nately, there are certain problems with this isotope. The 
first one is that the SU(3) symmetry is onl y ap proxi- 
mate as the scattering lengths are different |25|]. The 
more substantial problem, however, is that the three- 
body loss rate was found to be rather large in experi- 
ments [26j. This loss mechanism is not included in the 
model in Eq. ([2]). Although there are approaches where 
strong three-body loss is used as a dynamical three-body 
constraint [27| . a very important aspect, "baryon forma- 
tion" , cannot be studied that way. 

It is therefore desirable to find an isotope where three- 
body loss is much weaker. A possibility is to use 173 Yb, 
which was already trapped in an optical lattice (28j . Ac- 
tually, there is a simple argument why three-body re- 
combination should be less relevant than in 6 Li. In cases 
when the scattering length is much larger than the inter- 
atomic van der Waals length, the three-body recombina- 
tion rate K3 scales with some average scattering length 
a s and atomic mass to, 123. l30j 



(3) 
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As 6 Li is much lighter and has a larger scattering length 
than 173 Yb, the loss rates in 173 Yb are expected to be 
much weaker than in 6 Li. 

There is a further advantage of using the isotope 173 Yb: 
due to its closed electronic subshell, the state of the atom 
is determined by the nuclear-spin state. Since the inter- 
action between the atoms is mediated by dipolar van der 
Waals forces, it is in turn independent of the nuclear spin 
with a very good approximation. As a consequence, the 
interactions between 173 Yb atoms will have an increased 
symmetry, SU(2) — > SU(N), at ultracold temperatures 



31|. In Ref. [28j . SU(6) symmetry was established; how- 
ever, any number N < 6 of the nuclear states can be pop- 
ulated using optical pumping. The number N is fixed by 
the initial preparation due to the conservation of atoms 
in each component. From now on we shall concentrate 
on the case N = 3. 

The properties of 1 Yb make it a promising candi- 
date to realize the SU(3) Hubbard model: however, the 
scattering length is positive, a s ~ lOnm, [28( leading to 
repulsive interactions as U ~ a s . This is precisely a situ- 
ation when negative absolute temperatures would allow 
one to realize a system with effectively attracting inter- 
actions. 

Now we review some properties of the color supcrfluid 
and trionic phases based on Ref. [22| , where the half-filled 
SU(3) Hubbard model was studied at finite temperature. 
The entropy of the color superfluid can be seen in Fig. 2 of 
Ref. [22j ■ To enter this phase one requires an entropy per 
site S/L 3 rj 0.3 corresponding to an entropy per particle 
(with 3/2 atoms per site at half band-filling) s sa 0.2, at a 
temperature T w 0.2W, where W corresponds to the full 
bandwidth of the non-interacting band. Although the 
calculations were performed on a Bethc lattice with in- 



finite coordination number, z —¥ 00, the non-interacting 
density of states approximates the density of states for 
nearest-neighbor hopping in a simple cubic lattice. Thus 
the values listed here should approximately be valid for 
a d = 3 fcrmionic cloud in a cubic optical lattice. Re- 
garding the trionic phase, an artifact of the analysis on 
the Bcthe lattice with z — >• 00 is that the trions are lo- 
calized as their hopping scales with z~ 3 l 2 . Nevertheless, 
trions in d — 3 dimensions should compose some Fermi 
liquid which crosses over to the high-temperature phase 
of a three-component Fermi liquid. 

Unfortunately, the lowest value of the entropy per par- 
ticle reported for (two-component) fermions is so ~ In 2 
in the center of a trap (with Stot/N to t ~ 1.3 for the to- 
tal cloud) [H, (33)]. The main reason for these high en- 
tropics, i.e., hot fcrmionic clouds, is that Pauli blocking 
suppresses cooling efficiency at ultracold temperatures. 
New cooling procedures are being developed to reach sub- 
stantially lower entropies, e.g., to realize the Neel antifer- 
romagnetic state in the two-component Hubbard model 
[33l ] . In the following, however, we shall consider a "hot" 
fcrmionic cloud and expand the free energy in the pa- 
rameter ptF [HI . A reasonable agreement between high- 
temperature expansions and dynamical mean-field the- 
ory has been found in Ref. [32| for typical experimental 
parameters. 

We consider the free energy F = F(J3, fi, U, W) of the 
fermions, defined by 

- /3F = lnTTCnp[-f3(H su(3) - fiN + W^2u)], (4) 

i 

for a given chemical potential \i. Here we included also a 
three-body term ~ W to be able to determine the density 
of triple occupancies ti = fii\fii%hi3. We find that for a 
simple cubic lattice, the free-energy densit y / = F/L 3 
can be expanded to second order in pt F as [34| 



x 2 + <D(pt F ) 4 



(5) 



with 



zo = z {/3, ft, U, W) = 1 + 3C + X u + C u 3 w (6) 



and 

X 2 = X 2 (/3,fji,U,W) = 3C + 12C 2 

+6C 3 u MM 



1 - u 



1 — u 2 w 
(3(2U + W) 

with f = e^, u = e~@ u , and w = e~@ w . The relevant lo- 
cal quantities are calculated analytically in local-density 
approximation as derivatives of the free energy accord- 
ing to the usual definitions, n = —df/dfi, d = df/dU, 
t = df/dW, k = -df/dtp, and s = -df/dT, taken 
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at n = /j,o — VqT 2 and W = 0. These correspond to 
the local densities of total particle number, total double 
occupancy, triple occupancy, hopping amplitudes, and 
entropy, respectively. 

Now we turn to the protocol shown in Fig. [T] to reach 
the state with effectively attracting interactions at T < 0. 
The high-temperature expansion describes both the ini- 
tial and the final states of the cloud in an optical lattice, 
since it can be applied for both f3 > and j3 < 0. Ini- 
tially the system is in equilibrium with parameter val- 
ues Vq = O.Olip, U = t F , and f3 = 0.5/ip. We varied 
the central chemical potential [io so the total particle 
number N to t is changed. For the parameters being used, 
N to t ~ 5 x 10 4 corresponds to a metallic cloud with an 
entropy Stat /Ntot ~ 2.4, while N tot w 6.5 x 10 5 corre- 
sponds to a band insulator with Stot/N to t ~ 0.75. Thus 
realistic values of the entropy per particle are covered by 
the calculations. 

The final state corresponding to /?/ < is determined 
by energy (and particle) conservation. We approximate 
the density distributions at step 4 with the initial equi- 
librium distributions. Therefore we have to solve 



d 3 r [-t F (l - q)h{r) + Udi{v) - V r 2 ni {r)] 
d 3 r [-t F k f (r) + Ud f (r) - V r 2 n f (r)] ; 

d s rn i {r) = /rf 3 rn/(r), (8) 



U = t F 



for the final inverse temperature, /?/, and the final cen- 
tral chemical potential, [iQf. Here we used abbrevi- 
ations for initial and final quantities, e.g., the parti- 
cle density is nj(r) = n(r; (3, fi , Vq, U) and ri/(r) = 
n(r; Pf, [Xof, — Vq, U), respectively. The quantity q cor- 
responds to the relative dephasing of nearest-neighbor 
sites during step 3 [i7| : while the external potential is 
reversed, the cloud is held for some waiting time t w 
in a deepened optical lattice with a quenched hopping 
t F — > 0. During this step, each site acquires a differ- 
ent local phase due to the inhomogeneity induced by the 
external potential. When the optical lattice is finally re- 
laxed during step 4, the kinetic energy represents a sum 
of terms with random complex phases, which averages to 
zero if t w is large enough. In our calculations, we used 
two cases of no (q = 0) or complete (q = 100%) dephas- 
ing. 

After obtaining the parameters /3/ and ^o/ numeri- 
cally, we can calculate all quantities in the final state. 
We show the relative entropy change, Sf/Si — 1, the 
final inverse temperature |/?/|, and the relative "trion 
content" T to t / (N to t / 3) as a function of the compres- 
sion |Vb| (^r) 2 ^ 3 011 Fig. [3] (For an inhomogeneous 
system in a parabolic potential, the compression de- 
scribes the system in the thermodynamic limit, where 
N tot -+ oo,V -> 0, while V N^ t d is fixed @.) We see 
that deep in the band insulating regime, the entropy in- 
crease is around 30 % (or 60%) in the case of complete 



1 009r depilating 
1 P©- -O center, 0% d. 
> center, 100% 



5 10 
IV I (N tot l4K) m lt F 



Q-Q 0% dephasing 

^— O 100% dephasing 
qI — i — , — . — , — , — , — , — I — i — I 
3 6 9 12 15 
IV I (N tot /4n) m h F 



5« 0.5 - 
0.25 



O-0 0% dephasing 
0—0 100% dephasing _ 

- o o o o — e — e — e — e — e — e — o 



^- o ooo — o — o — o — o — o — o — o — o 



3 6 9 12 15 

IV I (N, 01 /47l) 2/3 /t F 

FIG. 3: (Color online) Lower panel: The relative change in 
the entropy as a function of the compression, after reversing 
Vo — > —Vo in the SU(3) Hubbard model for U = tF- Upper 
right panel: the final inverse temperature \/3f\. Upper left 
panel: Relative total trion occupation, T to t/{Ntot/3) in the 
final state is shown by the large symbols, while the relative 
trion occupation in the center, i(r = 0)/(n(r = 0)/3) is shown 
by the small symbols. 



(or no) dephasing. These values are not too high in com- 
parison to an ftH00% increase in the entropy of a two- 
component Fermi gas when it is first loaded in an op- 
tical lattice and then unloaded from it [32j . The final 
values of the entropy with the high trionic occupations 
should allow us to explore the crossover from the high- 
temperature phase to the trionic phase. In order to reach 
the color superfluid phase, however, additional reduction 
of the entropy is necessary. 



II. S=l ANTIFERROMAGNETIC SPIN CHAIN 

In our second example, we consider spin-1 bosons in a 
d = 1 dimensional lattice. Two ultracold bosonic atoms 
with spins S±, S2 = 1 can scatter in total spin sectors 
Si + S2 = and Si + S2 = 2. As a consequence, the 
appropriate Hamiltonian in a deep enough optical lattice 
in a homogeneous situation is [35l |36| 



(9) 



^^n,(n^l) + ^^(S l 2 -2n J ;), 



where ts > is the hopping amplitude; bi a are bosonic 
operators with spin projections a = +1,0,-1; hi = 
^2 a b\ a bia is the number of bosons at a site; Uq ~ 
(ao + 2a2) and U2 ~ (02 — o-o) ar e the interaction pa- 
rameters with as being the s-wave scattering lengths for 
scattering in the S = 0, 2 spin channels; and Si being 
the spin operator at site i. Typically, ao « 02, therefore 
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FIG. 4: (Color online) Schematic phase diagram of the d — 
1 quadratic-biquadratic spin model, after Ref. [13]. Arrows 
indicate parameters that apply to realizations with S7 Rb and 
23 Na isotopes. 



| U2 1 <C \Uq\. For ts <C Uq, a Mott insulating phase is re- 
alized when the filling is one boson per lattice site. At low 
energies in the Mott phase, one can describe the system 
by the so-called quadratic-biquadratic spin model [36H38j 



(ij) 



S,-' 



(10) 



The couplings are given by[3£ 



J1 



' B 
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3 U Q + U 2 3U - 2U 2 



(11) 



The model in Eq. (|10|) in d = 1 dimension has a rich 
ground state phase diagram, shown in Fig. @] (see also 
Ref. [37| and references therein; a phase diagram in a 
magnetic field can be found in Ref. |38[). In addition to 
these intriguing phases, even finite temperature (T =/= 0) 
properties are interesting. In the last decades, major 
theoretical efforts have been made to calculate T 7^ 
dynamical correlation functions of (gapped antifcrromag- 
netic) spin chains [39j . A simulation with ultracold atoms 
could provide an important experimental benchmark. 
However, since bosons collapse for Uq < 0, so far only the 
Ji, J 2 > regime could be realized. Negative absolute 
temperatures provide a possibility to reach the J\, J2 < 
regime. 

To realize the antiferromagnctic couplings, we revisit 
the approach of Ref. (3tJ and modify it to fit to the gen- 
eral protocol shown in Fig. [T] In Ref. [1?| constant elec- 
tric and staggered magnetic fields have been proposed to 
change the couplings and to reach the ground state of H 
from the ground state of — H . An energy gap is essential 
for the adiabatic evolution; i.e., the critical phase cannot 
be accessed. Application of constant electric fields is also 
not ideal for two reasons. First, there is no equilibrium 
for p 7^ in the thermodynamic limit in the presence 
of a constant electric field for a homogeneous system. 



Second, as discussed in Ref. [40(, a finite cloud in the 
presence of constant forces may exhibit a behavior that 
is substantially different from the corresponding homo- 
geneous system. Such dynamics should be avoided when 
the aim is to simulate a specific quantum system. 

Although we think that application of electric fields 
should be avoided, staggered magnetic fields are neces- 
sary for optimal results. To realize this staggered field, 
one can apply a secondary, spin-dependent optical lattice 
with a lattice period twice as large as the original which 
shifts the energies of the S z = ±1 atoms in opposite di- 
rections [13]. It can be represented by the term 



H st = -B^(-1)*SJ 

i 

= -b (-iMlA 



(12) 



i,cr = +l,0,-l 

The improved protocol to reach T < with spin-1 bosons 
involves taking both Vq —Vq and B — > —B during step 
3. This latter could be implemented experimentally by 
shifting the phase of the lasers that create the secondary 
optical lattice by ir. Note that in contrast to the case 
of spinless bosons, discussed in Ref. [13], the interaction 
Uo > is not to be changed. This repulsive interaction 
translates to effectively attracting bosons at < 0, and 
to antiferromagnctic couplings Ji, J 2 < in Eqns. ([TTJ) . 

We would like to emphasize that, in contrast to the 
instability of a spinor Bose gas prepared with as < 0, 
suddenly changing from the Mott insulator Uo ^> tg 
phase to the effectively attractive case Uq = — Uq < is 
metastable. Since only the atoms themselves can trans- 
port energy, the binding energy of two bosons has to be 
taken away by a large number w Uq/Ib of other atoms. 
Nevertheless, such a configuration has a very low prob- 
ability. As a consequence, double occupancies can form 
only at an exponentially suppressed rate. A similar dis- 
cussion was applied for a complementary situation in 
Ref. [U. 

Relaxation to an equilibrium state in a d = 1 dimen- 
sional system is a delicate and extremely difficult prob- 
lem, especially in the vicinity of points where the model 
is integrable |42j] . In fact, this is an additional reason why 
the experimental simulation of the spin model in Eq. (|T0|) 
is important. To resolve the problem of one-dimensional 
thermalization and to simulate the conditions in experi- 
ments, we consider an array of chains instead of a single 
chain: The bosons are most likely trapped initially in an 
anisotropic d = 3 dimensional optical lattice, where hop- 
ping along the chains, characterized by ts, is stronger 
than hopping between the chains, given by tg. Even 
weak interchain couplings should break integrability and 
allow relaxation to a thermal state. Note that with 
orthogonally arranged laser beams, ts and tg can be 
changed independently. Although one could also tune 
the aspect ratio of the parabolic potential, we focus on a 
spherically symmetric external potential. 

For the remainder of this section we focus on a cloud 
with interaction strengths Uq = 50ts, U2 = — O.OIL/q, but 
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similar results can be found for U2 = O.OAUn. These val- 
ues apply for 87 Rb and 23 Na, respectively [43 • Initially, 
the external potential is set to Vq = 0.005<s, and the cen- 
tral chemical potential fj-o is changed to vary the number 
of particles N to t and, consequently, the compression. The 
value of the interchain hopping is t B = 0At B . 

To describe the initial system, we use a Gutzwillcr 
wave function of the form 



i g > = n 



fata (^) 

~7T 



(air 



|0) 



(13) 



where f\(i) correspond to occupation amplitudes of local 
configurations, and the boson operators oj^ are related 
to the operators b\ a in Eq. ((9J by a canonical transfor- 
mation, 1361 



V2 



-,4 



V2 



(14) 



In Eq. (|13p we neglected triple- and higher occupancies 
which are almost completely suppressed in the strongly 
repulsing regime with at most one boson per lattice 
site. The Gutz wilier wave function works in the Mott 
insulating region well, while it is expected to shift the 
boundary of the Mott region in comparison to the ex- 
act value [l||. Since we focus on a tightly compressed 
cloud (no ~ Uo ^> t B ) in a relatively strong staggered 
field B > Ji, J2, the shifts of values of the total energies 
and total atom numbers (the quantities relevant for our 
calculations) are expected to be small. 

The variational energy per lattice site can be expressed 
as (with (G\G) = 1) 

E v = (G\H S1BH + H st - /i^n M |G)/L 3 

iot 

- t B Y^[Q* a (A)Q a (B) + Q* a (B)Q a (A)\ 

Ot 

-2t£ Y,^(A)Q a (B) + Q* a (B)Q a (A)} 

Ot 

En a (A) + n a (B) m A - m B 
9 B 9 



U d A + d B U 2 d A + d B 



(15) 



+ V2f: y (i)[f xx (i)-f yy (i)]-R. c .] 

dt = 2j2\U a (l)\ 2 + 2j2\Ml)\ 2 

a (3<a 

dl = -2E/-(0/M0 + 2El/^(0! 2 > (16) 

where I = A, B denotes sites in the two sublattices. One 
can reduce the number of variational parameters by using 
the sublattice symmetry. 

For given values of the control parameters one can min- 
imize Eq. (fTS")) with respect to the variational parameters 
f\. It is confirmed that double occupancies are negli- 
gible for the large value Uq = 50t B . The initial en- 
ergy is calculated using the densities obtained by the 
Gutzwiller ansatz in local-density approximation, e.g., 
n a (r) = n a (ji = Mo - V r 2 ), This gives 

E i (Uo,U 2 ,Vof = -V ,B f = -B) (17) 
= Jd 3 r |-t B (l + 2«)(1 - q)Y^[Q Aa {v)Q Ba {v) + h.c] 

(r) + n Ba (r)}/2 

Ot 

-B f [m A (r)-m B (r)]/2 
+ ^[d A (r)+d B (r)} + ^.[d A (r) + d B (v)\ 

Here k = t B /t B =0.1. Similarly to the fcrmionic case, 
the coefficient q corresponds to the dephasing efficiency 
during step 3. 

Let us now consider the final state. We assume that 
a high-temperature expansion [34l| in (3t B can be applied. 
The validity of this step has to be confirmed later. As we 
have reasoned above, all processes where a doubly occu- 
pied site is created are exponentially suppressed and will 
not occur during the typical time scale of an experiment. 
We therefore assume that the system relaxes in a sub- 
space where two-body occupations are absent. Therefore 
we use the following expansion of the free energy: 
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FIG. 5: (Color online) Final parameters for the S=l bosonic 
Hubbard model in a staggered field at a compression of 
\V \{3N tot /A7v) 2/:i /t B = 12 for t£ = O.lta, U = 50t s and 
IJi — — O.Olt/o- The final entropies and final inverse tem- 
peratures are shown for the cases of q = 0% and q = 100% 
dephasing of the initial kinetic energy. The different lines cor- 
respond to results calculated in 2nd and 4th-order expansions 
in /3tB- 



and 
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(22) 



where we introduced £ = e^, B = cosh((3B), and B = 
smh(/3B)/ (f3B). The densities of different quantities in 
local-density approximation [//(r) = fi — V r 2 ] can be cal- 
culated analytically from definitions like m s = —df/dB, 
fen = —df/dt B , and fcj^ = —df/dtg. Finally, the energy 
in the final state is given by 



E f = 



fd 3 r [- t B k\\(fif,iM,f,Bf,T) (23) 

-2t^k±H3f,nof,Bf,T) 
-B f m s (/3 f , Mo/, -B/, r) + V of r 2 n((3, (io,B f , r)] , 



where Bf = —B and Vqj = —Vq. Total-energy conser- 
vation implies that Ef equals the initial energy Ei given 
by Eq. ([T7|. Considering also the conservation of parti- 
cles Nf = Ni, one can find numerically the final inverse 
temperature f3f < and the final central chemical po- 
tential nof, which determine the final state. In Fig. [5] we 
show the final temperatures and the final entropies as a 
function of the staggered magnetic field for a given com- 
pression |V |(3iV tot /47r) 2 / 3 /is = 12, while in Fig. © we 
show quantities for B/ts — 0.3 as a function of the com- 
pression. As a check of the validity of the approach, we 
compare the results using the second- and fourth-order 
expansions. 



FIG. 6: (Color online) The final entropy of spin-1 bosons as 
a function of the compression for a staggered field value B = 
0.3tB with 100% dephasing of the initial kinetic energy, using 
fourth order expansion. Inset: the final inverse temperature 
13 f. Values of t B , Uo and Ui are the same as in Fig. [5] 



As we can see, for low values of B, and especially with 
no dephasing (q = 0) during step 3, the entropy per par- 
ticle is large, Stot/N to t > hi 3. This characteristic value 
corresponds to a system where sites arc populated by 
precisely one of three kind of particles in an uncorrclatcd 
fashion. Higher entropies indicate considerable fluctua- 
tions in not just the spin, but also in the occupation of 
the sites. This implies that the spin Hamiltonian defined 
in Eq. (fTUj) should not be used to describe the model, 
and that the application of the staggered magnetic field 
is necessary. As the strength of the staggered field is in- 
creased, both the entropy and \Tf\ decrease. However, 
the results of second- and fourth-order expansions start 
to deviate significantly for B/ts > 0.35. Further increas- 
ing B eventually leads to the breakdown of the high- 
temperature expansion. Although we cannot give quan- 
titative estimates based on our method, stronger stag- 
gered magnetic fields should lead to even lower entropies 
based on a simple physical picture. As H s t becomes the 
dominant term in the effective Hamiltonian, the initial 
densities at /3 — > +oo and B > approximate the densi- 
ties in the final state, with /3f < and Bf = —B, very 
well in the Mott insulating center. The heating occurs 
mainly at the compressible edges. It is important, how- 
ever, that B <§; Uo is maintained for metastability as the 
local Zccman energy of two atoms lias to be much smaller 
than the interaction energy. After step 5, Bf and t B have 
to be weakened adiabatically to reach the Hamiltonian in 
Eq. (|10p with antiferromagnetic couplings. We empha- 
size that exploring finite temperature dynamics of spin 
chains experimentally would be as interesting as estab- 
lishing the ground state phase diagram. 



III. CONCLUSIONS 

We have discussed two important theoretical models 
which can be simulated with an ultracold atomic cloud 
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at negative absolute temperatures. In the first case, re- 
pulsively interacting fermionic atoms are considered to 
realize the attractive SU(3) Hubbard model, while in the 
second case, bosonic atoms with spin S = 1 in the Mott 
phase at T < could simulate antiferromagnetically cou- 
pled spin chains. In general, considering negative tem- 
peratures provides an alternative way to realize couplings 
which are hard to reach at T > 0. 

We would like to emphasize that the models discussed 
here could be simulated in other ways. For example, in 
the case of the bosons one could reverse the interaction 
Uq — > — Uq using Fcshbach resonances during step 3 in- 
stead of reversing the external potential Vo- 
lt would be interesting to study the real-time dynam- 
ics of the relaxation in the two models. In the case of 



the fermionic cloud, we expect that energy diffusion de- 
termines the time scales, similarly to the two-component 
case [13, E^. The relaxation in the bosonic case is, how- 
ever, different and much more challenging. 
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